Genetic introgression from commercial European pigs to the indigenous Chinese Lijiang breed and associated changes in phenotypes.

BACKGROUND: Gene flow is crucial for enhancing economic traits of livestock. In China, breeders have used hybridization strategies for decades to improve livestock performance. Here, we performed whole-genome sequencing of a native Chinese Lijiang pig (LJP) breed. By integrating previously published data, we explored the genetic structure and introgression of genetic components from commercial European pigs (EP) into the LJP, and examined the impact of this introgression on phenotypic traits. RESULTS: Our analysis revealed significant introgression of EP breeds into the LJP and other domestic pig breeds in China. Using a haplotype-based approach, we quantified introgression levels and compared EP to LJP and other Chinese domestic pigs. The results show that EP introgression is widely prevalent in Chinese domestic pigs, although there are significant differences between breeds. We propose that LJP could potentially act as a mediator for the transmission of EP haplotypes. We also examined the correlation between EP introgression and the number of thoracic vertebrae in LJP and identified VRTN and STUM as candidate genes for this trait. CONCLUSIONS: Our study provides evidence of introgressed European haplotypes in the LJP breed and describes the potential role of EP introgression on phenotypic changes of this indigenous breed.

Deciphering the differential expression patterns of yield-related negative regulators in hexaploid wheat cultivars and hybrids at different growth stages.

BACKGROUND: Hexaploid bread wheat underwent a series of polyploidization events through interspecific hybridizations that conferred adaptive plasticity and resulted in duplication and neofunctionalization of major agronomic genes. The genetic architecture of polyploid wheat not only confers adaptive plasticity but also offers huge genetic diversity. However, the contribution of different gene copies (homeologs) encoded from different subgenomes (A, B, D) at different growth stages remained unexplored. METHODS: In this study, hybrid of elite cultivars of wheat were developed via reciprocal crosses (cytoplasm swapping) and phenotypically evaluated. We assessed differential expression profiles of yield-related negative regulators in these cultivars and their F1 hybrids and identified various cis-regulatory signatures by employing bioinformatics tools. Furthermore, the preferential expression patterns of the syntenic triads encoded from A, B, and D subgenomes were assessed to decipher their functional redundancy at six different growth stages. RESULTS: Hybrid progenies showed better heterosis such as up to 17% increase in the average number of grains and up to 50% increase in average thousand grains weight as compared to mid-parents. Based on the expression profiling, our results indicated significant dynamic transcriptional expression patterns, portraying the different homeolog-dominance at the same stage in the different cultivars and their hybrids. Albeit belonging to same syntenic triads, a dynamic trend was observed in the regulatory signatures of these genes that might be influencing their expression profiles. CONCLUSION: These findings can substantially contribute and provide insights for the selective introduction of better cultivars into traditional and hybrid breeding programs which can be harnessed for the improvement of future wheat.

Some plasma biomarkers of residual feed intake in beef cattle remain consistent regardless of intake level.

This study investigated whether plasma biomarkers of residual feed intake (RFI), identified under ad libitum feeding conditions in beef cattle, remained consistent during feed restriction. Sixty Charolais crossbred young bulls were divided into two groups for a crossover study. Group A was initially fed ad libitum (first test) and then restricted (second test) on the same diet, while Group B experienced the opposite sequence. Blood samples were collected from the 12 most divergent RFI animals in each group at the end of the first test and again after the second test. 12 plasma variables consistently increased, while three consistently decreased during feed restriction (FDR < 0.05). Only two metabolites, α-aminoadipic acid for Group A and 5-aminovaleric acid for Group B, were associated with RFI independent of feed intake level (FDR < 0.05), demonstrating moderate-to-high repeatability across feeding levels (intraclass correlation coefficient ≥ 0.59). Notably, both metabolites belong to the same metabolic pathway: lysine degradation. These metabolites consistently correlated with RFI, irrespective of fluctuations in feed intake, indicating a connection to individual metabolic processes influencing feed efficiency. These findings suggest that a portion of RFI phenotypic variance is inherent to an individual's metabolic efficiency beyond variations in feed intake.

Hybridization mediated range expansion and climate change resilience in two keystone tree species of boreal forests.

Current global climate change is expected to affect biodiversity negatively at all scales leading to mass biodiversity loss. Many studies have shown that the distribution of allele frequencies across a species' range is often influenced by specific genetic loci associated with local environmental variables. This association reflects local adaptation and allele changes at those loci could thereby contribute to the evolutionary response to climate change. However, predicting how species will adapt to climate change from this type of data alone remains challenging. In the present study, we combined exome capture sequences and environmental niche reconstruction, to test multiple methods for assessing local adaptation and climate resilience in two widely distributed conifers, Norway spruce and Siberian spruce. Both species are keystone species of the boreal forest and share a vast hybrid zone. We show that local adaptation in conifers can be detected through allele frequency variation, population-level ecological preferences, and historical niche movement. Moreover, we integrated genetic and ecological information into genetic offset predictive models to show that hybridization plays a central role in expanding the niche breadth of the two conifer species and may help both species to cope better with future changing climates. This joint genetic and ecological analysis also identified spruce populations that are at risk under current climate change.

Synthetic eco-evolutionary dynamics in simple molecular environment.

The understanding of eco-evolutionary dynamics, and in particular the mechanism of coexistence of species, is still fragmentary and in need of test bench model systems. To this aim we developed a variant of SELEX in vitro selection to study the evolution of a population of ∼1015 single-strand DNA oligonucleotide 'individuals'. We begin with a seed of random sequences which we select via affinity capture from ∼1012 DNA oligomers of fixed sequence ('resources') over which they compete. At each cycle ('generation'), the ecosystem is replenished via PCR amplification of survivors. Massive parallel sequencing indicates that across generations the variety of sequences ('species') drastically decreases, while some of them become populous and dominate the ecosystem. The simplicity of our approach, in which survival is granted by hybridization, enables a quantitative investigation of fitness through a statistical analysis of binding energies. We find that the strength of individual resource binding dominates the selection in the first generations, while inter- and intra-individual interactions become important in later stages, in parallel with the emergence of prototypical forms of mutualism and parasitism.

Ancient and recent origins of shared polymorphisms in yeast.

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

Portability of genomic predictions trained on sparse factorial designs across two maize silage breeding cycles.

KEY MESSAGE: We validated the efficiency of genomic predictions calibrated on sparse factorial training sets to predict the next generation of hybrids and tested different strategies for updating predictions along generations. Genomic selection offers new prospects for revisiting hybrid breeding schemes by replacing extensive phenotyping of individuals with genomic predictions. Finding the ideal design for training genomic prediction models is still an open question. Previous studies have shown promising predictive abilities using sparse factorial instead of tester-based training sets to predict single-cross hybrids from the same generation. This study aims to further investigate the use of factorials and their optimization to predict line general combining abilities (GCAs) and hybrid values across breeding cycles. It relies on two breeding cycles of a maize reciprocal genomic selection scheme involving multiparental connected reciprocal populations from flint and dent complementary heterotic groups selected for silage performances. Selection based on genomic predictions trained on a factorial design resulted in a significant genetic gain for dry matter yield in the new generation. Results confirmed the efficiency of sparse factorial training sets to predict candidate line GCAs and hybrid values across breeding cycles. Compared to a previous study based on the first generation, the advantage of factorial over tester training sets appeared lower across generations. Updating factorial training sets by adding single-cross hybrids between selected lines from the previous generation or a random subset of hybrids from the new generation both improved predictive abilities. The CDmean criterion helped determine the set of single-crosses to phenotype to update the training set efficiently. Our results validated the efficiency of sparse factorial designs for calibrating hybrid genomic prediction experimentally and showed the benefit of updating it along generations.

Integrating Dynamic 3D Chromatin Architecture and Gene Expression Alterations Reveal Heterosis in Brassica rapa.

Heterosis plays a significant role in enhancing variety, boosting yield, and raising economic value in crops, but the molecular mechanism is still unclear. We analyzed the transcriptomes and 3D genomes of a hybrid (F1) and its parents (w30 and 082). The analysis of the expression revealed a total of 485 specially expressed genes (SEGs), 173 differentially expressed genes (DEGs) above the parental expression level, more actively expressed genes, and up-regulated DEGs in the F1. Further study revealed that the DEGs detected in the F1 and its parents were mainly involved in the response to auxin, plant hormone signal transduction, DNA metabolic process, purine metabolism, starch, and sucrose metabolism, which suggested that these biological processes may play a crucial role in the heterosis of Brassica rapa. The analysis of 3D genome data revealed that hybrid F1 plants tend to contain more transcriptionally active A chromatin compartments after hybridization. Supplementaryly, the F1 had a smaller TAD (topologically associated domain) genome length, but the number was the highest, and the expression change in activated TAD was higher than that of repressed TAD. More specific TAD boundaries were detected between the parents and F1. Subsequently, 140 DEGs with genomic structural variants were selected as potential candidate genes. We found two DEGs with consistent expression changes in A/B compartments and TADs. Our findings suggested that genomic structural variants, such as TADs and A/B chromatin compartments, may affect gene expression and contribute to heterosis in Brassica rapa. This study provides further insight into the molecular mechanism of heterosis in Brassica rapa.

Novel [1,3,4]Thiadiazole[3,2-a]pyrimidin-5-ones as Promising Biofilm Dispersal Agents against Relevant Gram-Positive and Gram-Negative Pathogens.

Biofilm-associated infections pose significant challenges in healthcare settings due to their resistance to conventional antimicrobial therapies. In the last decade, the marine environment has been a precious source of bioactive molecules, including numerous derivatives with antibiofilm activity. In this study, we reported the synthesis and the biological evaluation of a new series of twenty-two thiadiazopyrimidinone derivatives obtained by using a hybridization approach combining relevant chemical features of two important classes of marine compounds: nortopsentin analogues and Essramycin derivatives. The synthesized compounds were in vitro tested for their ability to inhibit biofilm formation and to disrupt mature biofilm in various bacterial strains. Among the tested compounds, derivative 8j exhibited remarkable dispersal activity against preformed biofilms of relevant Gram-positive and Gram-negative pathogens, as well as towards the fungus Candida albicans, showing BIC50 values ranging from 17 to 40 µg/mL. Furthermore, compound 8j was in vivo assayed for its toxicity and the anti-infective effect in a Galleria mellonella model. The results revealed a promising combination of anti-infective properties and a favorable toxicity profile for the treatment of severe chronic biofilm-mediated infections.

The extent of introgression between incipient Clarkia species is determined by temporal environmental variation and mating system.

Introgression is pervasive across the tree of life but varies across taxa, geography, and genomic regions. However, the factors modulating this variation and how they may be affected by global change are not well understood. Here, we used 200 genomes and a 15-y site-specific environmental dataset to investigate the effects of environmental variation and mating system divergence on the magnitude of introgression between a recently diverged outcrosser-selfer pair of annual plants in the genus Clarkia. These sister taxa diverged very recently and subsequently came into secondary sympatry where they form replicated contact zones. Consistent with observations of other outcrosser-selfer pairs, we found that introgression was asymmetric between taxa, with substantially more introgression from the selfer to the outcrosser. This asymmetry was caused by a bias in the direction of initial F1 hybrid formation and subsequent backcrossing. We also found extensive variation in the outcrosser's admixture proportion among contact zones, which was predicted nearly entirely by interannual variance in spring precipitation. Greater fluctuations in spring precipitation resulted in higher admixture proportions, likely mediated by the effects of spring precipitation on the expression of traits that determine premating reproductive isolation. Climate-driven hybridization dynamics may be particularly affected by global change, potentially reshaping species boundaries and adaptation to novel environments.

The temporal and genomic scale of selection following hybridization.

Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the evolutionary dynamics within hybrid populations that underlie these patterns have been lacking. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of ancestry variation at varying spatial genomic scales through time. Here, we develop methods based on the Discrete Wavelet Transform to study the genomic scale of local ancestry variation and its association with recombination rates and show that these methods capture temporal dynamics of drift and genome-wide selection after hybridization. We apply these methods to published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio) and to inferred Neanderthal introgression in modern humans. Across systems, upward of 20% of variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. Signatures of selection at fine genomic scales suggest selection over longer time scales; however, we suggest that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from contiguous segments of genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available and can help shed light on generalities of the genomic consequences of interspecific hybridization.

Adaptive introgression of a visual preference gene.

Visual preferences are important drivers of mate choice and sexual selection, but little is known of how they evolve at the genetic level. In this study, we took advantage of the diversity of bright warning patterns displayed by Heliconius butterflies, which are also used during mate choice. Combining behavioral, population genomic, and expression analyses, we show that two Heliconius species have evolved the same preferences for red patterns by exchanging genetic material through hybridization. Neural expression of regucalcin1 correlates with visual preference across populations, and disruption of regucalcin1 with CRISPR-Cas9 impairs courtship toward conspecific females, providing a direct link between gene and behavior. Our results support a role for hybridization during behavioral evolution and show how visually guided behaviors contributing to adaptation and speciation are encoded within the genome.

Introduction of a diverse genetic background of Pyrus into Malus through intergeneric hybridization.

Wide hybridizations across species and genera have been employed to enhance agriculturally important traits in crops. Within the tribe Maleae of the Rosaceae family, different genera and species exhibit several traits useful for increasing diversity and gene pool through hybridization. This study aimed to develop and characterize intergeneric hybrid individuals between Malus and Pyrus. Through seed germination, shoot multiplication, and rooting in vitro, acclimatized seedlings showing vegetative growth on their own roots were obtained from crosses of Malus × domestica pollinated by Pyrus communis, P. bretschneideri, and the Pyrus interspecific hybrid (P. communis × P. pyrifolia). Comparative analysis of leaf morphology, flow cytometry, and molecular genotyping confirmed the hybrid status of the individuals. Genome-wide genotyping revealed that all the hybrid individuals inherited genomic fragments symmetrically from the Malus and Pyrus parents. To the best of our knowledge, this is the first report on the development of intergeneric hybrid seedlings between Malus × domestica and P. bretschneideri. Furthermore, the Pyrus interspecific hybrid individual served as a bridge plant for introducing the genetic background of P. pyrifolia into Malus × domestica. The results of this study provided a crucial foundation for breeding through intergeneric hybridization between Malus and Pyrus, facilitating the incorporation of valuable traits from diverse gene pools.

Integrative taxonomic analyses reveal that rapid genetic divergence drives Acropora speciation.

Reef-building corals provide the structural basis for one of Earth's most spectacular and diverse but increasingly threatened ecosystems. The reef-building coral genus Acropora may have undergone substantial speciation during the Pleistocene climate and sea-level changes. Here, we aimed to evaluate the speciation history of four morphologically similar tabular Acropora species (Acropora aff. hyacinthus, A. cf. bifurcata, A. cf. cytherea, and A. cf. subulata) using an integrative approach with morphology, genetic, and reproduction methodology. Extensive morphological analyses showed that these four species are distinct and exhibited high gamete incompatibility, preventing hybridization. Furthermore, population structure and principal component analyses with SNPs (>60,000) indicated that these species were genetically distinct, and the ABBA-BABA test did not support introgression among these species. Many of their coding and noncoding RNA sequences showed high genetic variance at loci with high Fst values along the genome. Comparison of these orthologs with those of other Acropora species suggested that many of these genes are under positive selection, which could be associated with spawning time, gamete, and morphological divergence. Our findings show that the speciation of tabular Acropora occurred without hybridization, and the divergence accompanying the rapid evolution of genes in species-rich Acropora could be associated with speciation.

Emerging Trends in Wheat (Triticum spp.) Breeding: Implications for the Future.

Wheat (Triticum spp and, particularly, T. aestivum L.) is an essential cereal with increased human and animal nutritional demand. Therefore, there is a need to enhance wheat yield and genetic gain using modern breeding technologies alongside proven methods to achieve the necessary increases in productivity. These modern technologies will allow breeders to develop improved wheat cultivars more quickly and efficiently. This review aims to highlight the emerging technological trends used worldwide in wheat breeding, with a focus on enhancing wheat yield. The key technologies for introducing variation (hybridization among the species, synthetic wheat, and hybridization; genetically modified wheat; transgenic and gene-edited), inbreeding (double haploid (DH) and speed breeding (SB)), selection and evaluation (marker-assisted selection (MAS), genomic selection (GS), and machine learning (ML)) and hybrid wheat are discussed to highlight the current opportunities in wheat breeding and for the development of future wheat cultivars.

Application of machine learning for identification of heterotic groups in sunflower through combined approach of phenotyping, genotyping and protein profiling.

Application of machine learning in plant breeding is a recent concept, that has to be optimized for precise utilization in the breeding program of high yielding crop plants. Identification and efficient utilization of heterotic grouping pattern aided with machine learning approaches is of utmost importance in hybrid cultivar breeding as it can save time and resources required to breed a new plant hybrid/variety. In the present study, 109 genotypes of sunflower were investigated at morphological, biochemical (SDS-PAGE) and molecular levels (through micro-satellites (SSR) markers) for heterotic grouping. All the three datasets were combined, scaled, and subjected to unsupervised machine learning algorithms, i.e., Hierarchical clustering, K-means clustering and hybrid clustering algorithm (hierarchical + K-means) for assessment of efficiency and resolution power of these algorithms in practical plant breeding for heterotic grouping identification. Following the application of machine learning unsupervised clustering approach, two major groups were identified in the studied sunflower germplasm, and further classification revealed six smaller classes in each major group through hierarchical and hybrid clustering approach. Due to high resolution, obtained in hierarchical clustering, classification achieved through this algorithm was further used for selection of potential parents. One genotype from each smaller group was selected based on the maximum seed yield potential and hybridized in a line × tester mating design producing 36 F1 cross combinations. These F1s along with their parents were studied in open field conditions for validating the efficacy of identified heterotic groups in sunflowers genetic material under study. Data for 11 agronomic and qualitative traits were recorded. These 36 F1 combinations were tested for their combining ability (General/Specific), heterosis, genotypic and phenotypic correlation and path analysis. Results suggested that F1 hybrids performed better for all the traits under investigation than their respective parents. Findings of the study validated the use of machine learning approaches in practical plant breeding; however, more accurate and robust clustering algorithms need to be developed to handle the data noisiness of open field experiments.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression.

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

The emerging importance of cross-ploidy hybridisation and introgression.

Natural hybridisation is now recognised as pervasive in its occurrence across the Tree of Life. Resurgent interest in natural hybridisation fuelled by developments in genomics has led to an improved understanding of the genetic factors that promote or prevent species cross-mating. Despite this body of work overturning many widely held assumptions about the genetic barriers to hybridisation, it is still widely thought that ploidy differences between species will be an absolute barrier to hybridisation and introgression. Here, we revisit this assumption, reviewing findings from surveys of polyploidy and hybridisation in the wild. In a case study in the British flora, 203 hybrids representing 35% of hybrids with suitable data have formed via cross-ploidy matings, while a wider literature search revealed 59 studies (56 in plants and 3 in animals) in which cross-ploidy hybridisation has been confirmed with genetic data. These results show cross-ploidy hybridisation is readily overlooked, and potentially common in some groups. General findings from these studies include strong directionality of hybridisation, with introgression usually towards the higher ploidy parent, and cross-ploidy hybridisation being more likely to involve allopolyploids than autopolyploids. Evidence for adaptive introgression across a ploidy barrier and cases of cross-ploidy hybrid speciation shows the potential for important evolutionary outcomes.

Dysfunction of duplicated pair rice histone acetyltransferases causes segregation distortion and an interspecific reproductive barrier.

Postzygotic reproductive isolation, which results in the irreversible divergence of species, is commonly accompanied by hybrid sterility, necrosis/weakness, or lethality in the F1 or other offspring generations. Here we show that the loss of function of HWS1 and HWS2, a couple of duplicated paralogs, together confer complete interspecific incompatibility between Asian and African rice. Both of these non-Mendelian determinants encode the putative Esa1-associated factor 6 (EAF6) protein, which functions as a characteristic subunit of the histone H4 acetyltransferase complex regulating transcriptional activation via genome-wide histone modification. The proliferating tapetum and inappropriate polar nuclei arrangement cause defective pollen and seeds in F2 hybrid offspring due to the recombinant HWS1/2-mediated misregulation of vitamin (biotin and thiamine) metabolism and lipid synthesis. Evolutionary analysis of HWS1/2 suggests that this gene pair has undergone incomplete lineage sorting (ILS) and multiple gene duplication events during speciation. Our findings have not only uncovered a pair of speciation genes that control hybrid breakdown but also illustrate a passive mechanism that could be scaled up and used in the guidance and optimization of hybrid breeding applications for distant hybridization.

Genetic diversity based on multivariate techniques in elephant grass genotypes for bioenergy.

The use of biomass for energy production constitutes a promising strategy that warrants the search for new sources of biomass. Elephant grass has been gaining notoriety due to its high dry matter yield and rapid growth. The present study was carried out to quantify the genetic divergence of nine elephant grass half-sib families in order to identify genotypes with greater genetic divergence and productive potential for hybridization, using the hierarchical clustering methodology based on principal components. Half-sib families were generated using genotypes from the Active Germplasm Bank of Elephant Grass. The experiment was laid out in a randomized-block design with nine half-sib families, three replicates, and eight plants per plot. A total of 216 genotypes of elephant grass were evaluated. Principal component (PC), biplot, and hierarchical clustering analyses for diversity estimation were conducted using R software. The first two PCs of biplot analysis accounted for 64% of the cumulative variation. Dry matter yield was the most important trait for genotype discrimination (0.89), followed by plant height (0.67) and stem diameter (0.61) in PC1. In this analysis, the distances between accessions were considered and there were no family links, which indicates the existence of wide variability within the evaluated families, since genotypes belonging to the same family were not grouped together, but rather distributed into different groups. Crosses between genotypes of group three and genotypes of groups one and two are recommended for the development of high-yielding genotypes when aiming at energy production.